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Abstract 

Using the theory of elasticity of polymer gels we show that large- 
scale cross-link density patterns written into the structure of the net¬ 
work in the melt state, can be revealed upon swelling by monitoring 
the monomer density patterns. We find that while isotropic deforma¬ 
tions in good solvent yield magnified images of the original pattern, 
anisotropic deformations distort the image (both types of deforma¬ 
tion yield affinely stretched images in 9 solvents). We show that in 
ordinary solids with spatially inhomogeneous profile of the shear mod¬ 
ulus, isotropic stretching leads to distorted density image of this pro¬ 
file under isotropic deformation. Using simple physical arguments we 
demonstrate that the different response to isotropic stretching stems 
from fundamental differences between the theory of elasticity of solids 
and that of gels. Possible tests of our predictions and some potential 
applications are discussed. 
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1 Introduction 


Polymer networks are unique soft solids which can be significantly deformed 
without irreversible damage to their structure. A network is formed by cross- 
linking a melt or a semidilute solution of polymer chains. Once a homoge¬ 
neous (on length scales large compared to its “mesh” size) network is formed, 
one can generate large-scale patterns in it by further cross-linking, followed 
by swelling (and possibly stretching) of the network, resulting in a gel in- 
homogeously swollen by solvent. This can be done, for example, by adding 
light-sensitive cross-links to a transparent network. Focusing a laser beam in 
regions inside the gel one can “write” information into gel structure in the 
form of 2D or 3 D patterns of cross-linking density. In this paper we show 
that although such information is hidden at preparation conditions, it can 
be recovered by swelling the gel since unobservable variations of cross-link 
density in the melt are transformed into observable variations of monomer 
density in the swollen gel. 

Regions of a gel with increased cross-link concentration can be consid¬ 
ered as inclusions with enhanced elastic modulus. If such inclusions deform 
differently from polymer matrix, as in case of any normal elastic solids, they 
would induce elastic stresses in the gel and initial pattern would be signif¬ 
icantly distorted due to long range character of elastic interactions. This 
scenario determines, for example, the elastic properties of amorphous poly- 
crystalline solids but it does not apply to polymer gels, because of the unusual 
character of gel elasticity. We show that in swollen gels that are isotropically 
stretched by absorption of solvent, the observed monomer density pattern is 
not distorted and is simply an affinely stretched variant of the initial cross- 
linking pattern. Such gels can serve as a magnifying glass that enlarges the 
initially written pattern without distorting its shape. The corresponding 
magnification factor can be very large in case of super-elasic networks. 
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2 Free energy of a gel with cross-linking den¬ 
sity pattern 


In this paper we use the simplest mean field model of a gel with free energy 

m- 

A = J ^ G (*o) + f [p ( X( ))] j fix ° (!) 

Here p (xo) is monomer density as function of coordinates Xo in prepara¬ 
tion state. We assume that the gel was initially cross-linked in a polymer 
melt and then a pre-programmed pattern in cross-link concentration (i.e., 
a well-defined region of higher cross-link density compared to that of the 
surrounding network) is created in the network using, say, a light-sensitive 
cross-linking technique (the case of cross-linking in semi-dilute solution in 
good solvent is analyzed in SI). Here G (x 0 ) is the polymer contribution to 
the elastic modulus of the cross-linked melt (which is proportional to the 
local cross-link density) 

G (x 0 ) = G + G (x 0 ) (2) 


and G (x 0 ) represents the variations of cross-link density introduced by the 
second cross-linking step (Fig. HJi). / (p) is the osmotic (interaction) part of 
the free energy of the gel, with monomer density p. F is the deformation 
gradient tensor 


F = fai 
dx 


o j 


(3) 


and x are coordinates of deformed gel. It is convenient to assume that the 
gel is deformed with respect to preparation state in two stages: 


i 


(4) 


Thus, the gel is first stretched with respect to preparation state by factors \ 
along axes i. For such a deformation x = A • x 0 with components x* = A iX 0i 
and we get 

fy = A Aj, P = P = jrr (5) 

1L 'U 

where po is th e uniform monomer density in the undeformed state of prepa¬ 
ration (Fig. [T}}). Notice that the coordinates x describe a stretched network 
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Figure 1: Initial cross-link concentration c (x, y) a); monomer density profiles 
p(x,y ): of a gel in a melt state b), in reference state after stretching c), d), 
in stretched equilibrium state e),f). Gel is isotropically stretched by factors 
Ai = A 2 = 1.5 c), e) and anysotropically stretched by factors Ai = 2, A 2 = 1 


d),f). 


with inhomogeneous cross-link density but a homogeneous monomer density 
(Figs, dh and d). 

Even though such a homogeneous (in monomer density) state does not 
minimize the free energy and therefore is not an equilibrium state of the 
deformed gel, we use it as a reference state. The true equilibrium state of 
the deformed network has an inhomogeneous monomer density prohle and is 
defined by introducing a displacement field u (x) defined with respect to the 
above reference state: 


x.' = x + u (x) 


( 6 ) 


and we get gradient tensor and monomer density as function of coordinates 
x 



(7) 


( 8 ) 


Minimizing the free energy in Eq. (jT]) with respect to displacements u 
at the preparation state (all Aj = 1) we conclude that in a melt the cross¬ 
links and the monomers will remain at their previous possitions and the 
elastic reference state will not change after relaxation. We conclude that 
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information about the pattern written on network structure is hidden in 
preparation state and can only be revealed after swelling. 


3 What is the equilibrium density profile? 


In a swollen state the monomer density is small and the interaction energy 
can be expanded as / (p) ~ ksTBp 2 / 2, where k B is Boltzmann constant, T 
is temperature and B is second virial coefficient. Expanding the free energy 
in powers of u and integrating over the volume of the undeformed network 
with measure dx 0 = dx/ A* we get 

2 


AA = 




K„ 


dui 
^ dX; 


dxi 

dx 

TP? 


dxi 


K os = k B TBp 2 . 


(9) 


The equilibrium deformation of the gel is found by minimizing this free 
energy. Its variation is 


SA = 


Ji-W 


i d 2 Ui 

3 dx 2 


5ui 


ix ^ 2 Uj x„. 

K ° s / r\ r\ Oil} 

OXiOXj 
ij 1 


dx 

rii ^i 


and therefore, the minimum condition is 




'dx, 


dxidxn 


We are interested only in variations of monomer density 


pW 

p 


-E 


duj (x) 
dxi 


( 10 ) 


( 11 ) 


( 12 ) 


where p is average density. Taking the gradient of both sides of Eq. (HIT) we 
obtain an equation for the variations of monomer density 

d 2 p (x) 


(A 1 -x) 

ox 2 ^ 


7 i 


dx 2 


= 0 


P 


(13) 
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where 


( 14 ) 



The solution of this equation 


d( x ) 



/ 


g[ 7 1 • ( x - y)] G (A ^x) (15) 

1 I j fj 


is expressed through Green’s functions of the Laplace equation in 2 and 3 
dimensions, respectively: 



In case of isotropically stretched/swollen gel with all A* = A the equilib¬ 
rium monomer density depends on local cross-link concentration, 



(17) 


We conclude that under isotropic deformation such as swelling, the monomer 
density produces an undistorted, uniformly stretched image of the pattern of 
cross-link density originally “written” on the homogeneous network (compare 
Figs. [1^ and e). 

Equilibrium displacement is expressed through the variation of monomer 
density, Eq. (1171) . as 



(18) 


We conclude that although density variations in isotropically deformed gels 
are strictly local, there is long-range strain field decaying as power law of 
the distance |x — y|. This strain induces a stress distribution in the gel, 
which can be observed by measuring the birefringence of transmitted light 
(stress-optical law [3]). 

In anisotropically deformed networks the pattern is strongly distorted 
(compare Figs. [TJi and f) and p (x) decays as power law of a distance |x — y| 
from the localized cross-link density inhomogeneity G (y). Observe that vari¬ 
ations of monomer density are largest along the direction of stretching. This 
effect is closely related to the well known “butterfly” picture in contour plots 
of neutron scattering from random inhomogeneities of network structure in 
anisotropely deformed swollen gels [[Tj. 
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4 What is the difference between deforma¬ 
tion of gels and elastic solids? 


In order to understand the difference between gels and normal solids we recall 
that the free energy of any solid is a functional of the nonlinear strain tensor 
Uik 0 , 

dui . duk \ - dui duk 


2w ik I'ii l /.■ 


ij r kj 


dik 


dxk ^ dx 




dxj dxj 


(19) 


While the last term is usually neglected in the linear theory elasticity of 
solids because solids behave elastically only under small deformations, it can 
be shown that only this nonlinear term contributes to the elasticity of gels 
and that the elastic part of the free energy of gels (Eq. ([I])) is linear in this 
nonlinear strain W- Physically, the difference between elastic energy of a 
solid, which is a quadratic form in the linear strain and of a gel, which is linear 
in the nonlinear strain tensor, stems from the fact that while in solids there 
is a stress-free state of equilibrium (crystal lattice) that minimizes the energy 
of interaction between the atoms, the equilibrium state of gels is not stress- 
free. Polymer networks are made of entropic springs and, in the absence 
of osmotic pressure due to permeation by good solvent or due to excluded 
volume interactions in the melt state, such networks would collapse to the 
size of a single spring. The finite length of entropic springs in the swollen gel 
is the result of osmotic pressure which can be replaced by equivalent isotropic 
stretching forces that act on the outer boundaries of the gel |8]. 

The difference between gels and solids becomes apparent when consid¬ 
ering two simple toy models of heterogeneous gel and solid as two Hookean 
springs with moduli k\ and k 2 , connected in series as in Fig. [2j 

a) Gel model: Osmotic pressure is represented by a force f os applied to free 
ends of the connected springs. In the presence of this force the equilibrium 
lengths of the Gaussian springs become 7q q = f os /k i and r^ 1 = / os //c 2 , and 
if we apply additional force /, each of the springs will deform affinally with 
distance rq + r 2 between the ends to which the force is applied (’’boundaries” 
of the system): 


n = Arq q , r 2 = A r e 2 q , A = ^ ^ (20) 

rq + rq 

b) Solid model: The springs of a “solid” have equilibrium lengths r) q and 
rff 1 in the stress-free state. During stretching due to force / applied to the 
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Figure 2: Two springs models, demonstrating affine deformation of gels a) 
and non-affine deformation of solids b). 


ends of the two-spring system, such a solid deforms nonaffinelly: 

k 2 


r\ = r® q + (A — 1) (r® q + r. 2 q ) 
r 2 = r^ + (A-l)(r^ + r 2 eq ) 


k\ + k 2 
k\ 

k\ + k 2 


( 21 ) 


with the soft spring ( k\ < k 2 ) stretched more than the rigid one. 

These two simple toy models illustrate why under isotropic deformations, 
cross-linking density patterns in gels are stretched affinally, whereas soft re¬ 
gions in solids would undergo larger deformation compared to more rigid 
regions, thus distorting the original pattern. 


5 Discussion 

We studied the combined effect of swelling and deformation on inhomo¬ 
geneous networks, prepared by cross-linking a melt of polymer chains. It 
is well-known that cross-link density heterogeneities that have no effect on 
the monomer density in the state of preparation (a melt or a concentrated 
polymer solution), can be revealed by swelling the gel and observing the en¬ 
hancement of light, x-ray and neutron scattering from the resulting monomer 
density inhomogeneities pdl2]. In this paper we focused on a related phe¬ 
nomenon, namely that when large-scale cross-link density patterns are writ¬ 
ten into the network structure, the hidden image can be revealed by swelling 
and stretching the gel and observing the corresponding patterns of monomer 
density. Using the mean field theory of elasticity of polymer gels we showed 
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that stretching/swelling in good solvent acts as a magnifying glass: while 
isotropic stretching reproduces an enlarged but otherwise undistorted ver¬ 
sion of the original pattern, anisotropic stretching distorts this pattern, see 
figure m 

We compared these results with those obtained for ordinary elastic solids 
with inhomogeneous elastic moduli and found that in this case even isotropic 
deformations lead to distorted patterns. We showed that the fundamental 
difference between response of inhomogeneous gels and solids to isotropic 
stretching can be traced back to the fact that unlike regular springs that 
have an equilibrium length even in the absence of stress, the equilibrium 
length of entropic springs is entirely determined by the osmotic forces that 
isotropically stretch the polymer gel. 

Finally, we would like to comment on the possibility of experimental 
verification and on possible applications of our results. In most application 
involving gels such as biomimetic sensors, actuators and artificial muscles 
|13j . macroscopically inhomogeneous (layered) gels undergo shape transitions 
when the thermodynamic conditions are changed or in response to application 
of external fields p3]. In our case, the cross-link density pattern imprinted 
into the gel structure by, say, activation of light-sensitive cross-links, can be 
microscopic (micron size) and therefore would have little effect on the shape 
of the gel. Upon swelling and/or isotropic stretching in good solvent, the 
magnified density pattern can be imaged on a light-sensitive screen. The 
contrast can be significantly enhanced by stretching the entire gel in poor 
solvent or by focusing a laser beam on the localized pattern and heating it, 
resulting in local change of the quality of solvent. Finally, the sensitivity of 
the image to quality of solvent (the distortion under anisotropic deformation 
disappears in 9 solvent - see SI) can be useful for sensor devices. 
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Supporting Information. 

In the SI we show that patterns obtained by cross-linking a semi-dilute poly¬ 
mer solution, deform affinely (non-affinely) under isotropic (anisotropic) de¬ 
formation, just like in the case of cross-linking in the melt. We then ana¬ 
lyze how the pattern deforms under several different solvent conditions. We 
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show that the pattern always stretches affinely in a 0-solvent, even under 
anisotropic deformations. Since the contrast between the high and the low 
monomer density regions can be significantly enhanced in a poor solvent we 
proceed to analyze the density profiles in gels that are isotropically stretched 
in mildly poor solvents (at lower solubility, such stretched gels will undergo a 
transition into a strongly inhomogeneous state characterized by the appear¬ 
ance of dense filamentous structures m- We find that when the amplitude 
of cross-link density variations is sufficiently low, the image stretches affinely 
with the isotropic deformation but that for larger density contrasts the image 
becomes distorted, especially near the edges and corners of the pattern. 
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Supplementary Information to: 

Cross-linking patterns and their images in 
swollen and deformed gels 

SI 1. Gel free energy in a good solvent 

Consider a gel prepared in a good solvent at the monomer density p 0 that 
is swollen to density p. Its free energy is the sum of elastic and osmotic 
contributions. The osmotic pressure 7r of the gel in a good solvent increases 
proportionally to the 9/4 power of monomer density p [I] 

7T ~ ( k B T/b 3 ) (pfc 3 ) 9/4 (SI) 

where kg is Boltzmann constant, T is temperature and b is monomer size. 
The osmotic part A os of the free energy per polymer chain between network 
junctions is proportional to the free energy density (~ n) divided by the 
number of chains per unit volume ( p/N ), where N is the chain degree of 
polymerization: 

At- 7 ,/(p/N) ~k B TN(pb 3 ) 5/1 (S2) 

The dimension of the chain along the main axis a of deformation is R a = 
A a R 0 , where A Q is deformation factor along this axis (defined as eigenvalue 
of local deformation gradient tensor F) and i?o is the chain size in the state 
at which the gel was formed. The elastic free energy per chain is 



where Rfi is the amplitude of fluctuations of the chain in the deformed state. 
In a heterogeneous network the direction of the triad of deformation axes a 
depends on its position, and the sum of squares of local deformation factors 
in Eq. (lS3j) can be rewritten through the deformation gradient tensor F as 

E a « = E f «- < S4 ) 

a ij 

Since the mean-square amplitude of chain fluctuations is proportional to the 
mean-square polymer size at semi-dilute good solvent conditions and scales 
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with monomer density as |TJ R'jj ~ p 4 / 4 , while the mean-square chain size 

in the preparation conditions scales as ~ p q 1 ^ 4 , the elastic free energy per 
chain is 



(S5) 


At the equilibrium swelling (p = p eq ) in the absence of additional defor¬ 
mations the total free energy of the gel per chain is: 

N(p eq b 3 ) 5/4 + (po/p eq ) 5/12 j (S6) 


A dl = A% + A%~ k B T 


and it is minimized at the density [2] 


Peq ^ 


(Po& : 


3\l/4 


6 3 AT3/5 

corresponding to maximum swelling ratio 

A eq =(po/p eq ) 1 / 3 ^(p 0 fc 3 ) 1 / 4 iV 1/5 


(57) 

(58) 


Note that similar expression for \ eq is obtained in mean held model of a 
gel with second virial coefficient B ~ b 3 , see main text. This conclusion can 
also be extended to our solution of the image storing problem. Since both 
elastic (Eq. (jS5]) ) and osmotic (Eq. (jS2ll ) terms in the gel free energy are 
multiplied by the same scaling factor (p/po) 1 ^ 4 such scaling renormalization 
does not change the results obtained for the mean held model. 


SI 2. How does the pattern change in different 
solvent conditions? 

^-solvent 

In a 0-solvent the second virial coefficient vanishes (B = 0) and equation 

^ = Mx-y)]G(A-‘. x)l g- (SB, 

reproduces without distortion affinely stretched initial pattern 

P( x ) /P = G (A -1 - x ) /G, (S10) 

even for anisotropically stretched gels (small deviations from affinity are ex¬ 
pected because of the non-vanishing third virial coefficient) 
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Poor solvent 


Strong enhancement of monomer density contrast can be obtained by placing 
the gel (with fixed boundaries - otherwise it would collapse) in a poor solvent 
with negative second virial coefficient B < 0. In case of very poor solvent 
with 

7 2 = GX 2 + k B TBp 2 < 0 (Sll) 

the gel becomes unstable with respect to formation of domains with different 
monomer density [3]. Below we consider the case of poor solvent close to 
^-conditions with small B < 0 and positive 7 2 > 0. 

At small 7 2 the amplitude of density variations p (x) can be significantly 
increased because of the small denominator in equation 

and we have to take into account corrections due to second order in u term 



dui duj 
dxi dxj 


dui diij 
dxj dxi 


(S13) 


in expression for monomer density, 

P(x) = 

where 


P 


det (i^) ’ 


dui 


<M*3) = i + Eay + A < u 


To first order in A we find 

p(x) _ A 2 r 


P 
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G ( - ) - GA (u) 


(514) 

(515) 


(S16) 


where the equilibrium displacement in A (u) is determined as 

^( x )--^^-/0( x - y )G(|)dy (S17) 

We conclude that the correction term in Eq. (IS 161) enhances the contrast 
between the high and the low monomer density regions of the profile f Fig. IS 1 [) . 
This effect is the mostly pronounced near the corners of the pattern where 
several edges converge and it leads to distortion of the otherwise affinely 
stretched profile at these points. 
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Figure SI: Density profile in poor solvent for initial cross-link concentration 
shown in Fig. la). 


SI 3. Random heterogeneities 

Frozen-in random heterogeneities of network structure can change the image 
beyond recognition [4], The free energy of a gel with frozen-in heterogeities 
was derived in Ref. [5]. The only source of heterogeneities in the melt with 
fixed monomer density is statistical distribution of cross-links in the state of 
preparation that arises as the consequence of the random process of cross- 
linking. This frozen-in distribution is described by an additional contribution 
to the free energy: 



where / y - (x) is random Gaussian function of coordinate x, characterized by 
correlation function 

fij (x) hi (x') ~ ( k B T ) 2 05 (x - x') 5 ik 5 jt (S19) 

Comparing the amplitude of frozen-in fluctuations on a scale R with varia¬ 
tions of elastic modulus G on this scale we conclude that the contribution of 
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frozen-in heterogeneities can be neglected if 


G/G> 1/ (ci? 3 ) 1/2 (S20) 

where c is average cross-link concentration and thus, frozen-in heterogeities 
have no influence on large-scale patterns. The suppression of frozen-in het¬ 
erogeneities of monomer density is due to strong overlap of network chains [Bj . 
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